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Effective field theory methods are used to study factorization of the deep inelastic scattering cross- 
section. The cross-section is shown to factor in QCD, even though it does not factor in perturbation 
theory for some choices of the infrared regulator. Messenger modes are not required in soft-collinear 
effective theory for deep inelastic scattering as x —* 1. 


Many high energy processes in QCD factor into the 
convolution of a perturbatively calculable short dis¬ 
tance contribution,and a non-perturbative long distance 
hadronic function |l|. The classic example is deep inelas¬ 
tic scattering (DIS), where the structure functions fac¬ 
tor into short-distance hard-scattering kernels which de¬ 
pend on the large momentum scale Q, and parton distri¬ 
bution functions which depend on the non-perturbative 
scale Aqcd- DIS as x —> 1 has been studied using 
soft-collinear effective theory (SCET OBDD , which 
allows one to systematically separate out the different 
momentum scales, and also sum the Sudakov double¬ 
logarithms H0 

Recently, it was argued Jj that DIS does not factor 
as x —> 1 because of messenger modes jjj) of momentum 
~ Aqcd \J Aqcd / Q, which cause the hadronic function 
to depend on Q as well as Aqcd • The analysis of Ref. 0 
is based on a study of momentum regions contributing to 
perturbation theory graphs JTfj- The same method has 
been used in other examples, such as the calculation of 
the B —> 7r form factor, where it was also claimed that 
messenger modes are important 0 

The above claims are surprising. One would not expect 
the scale Aqcd \J Aqcd /Q to be relevant for hadronic 
structure functions and form factors. DIS has been exten¬ 
sively studied by diagrammatic methods, and has been 
shown to factor as x —> 1 into the convolution of a hard 
kernel, jet function, and parton distribution function, 
which depend on the scales Q , Q\f 1 — x and Aqcd, re¬ 
spectively [l]. The conclusions of Ref. *8] contradict these 
results. 

This paper discusses factorization and the (non) ex¬ 
istence of messenger modes for DIS. It is shown that 
a naive evaluation of diagrams in perturbation theory 
with off-shell external states leads to the conclusion that 
messenger modes exist, and that structure functions do 
not factor, as found in Ref. Jsj. However, a confining 
gauge theory such as QCD does not require messenger 
modes, and factorization holds. The basic flaw in the 
argument is equating momentum regions that contribute 
to a Feynman diagram in perturbation theory with prop¬ 
agating modes in an effective field theory description of 
a confining gauge theory such as QCD. Infrared modes 
which exist in perturbation theory are not present due 
to confinement. In an effective theory such as SCET, 
perturbation theory is used to sum logarithms using the 


renormalization group, and to compute short distance 
matching conditions. The infrared part of both the full 
and effective theory amplitudes is not given by pertur¬ 
bation theory, but nevertheless, perturbative amplitudes 
can be used to compute the matchin g co ndition, which 
only depend on short distance physics [Ty] • Perturbation 
theory is not used to compute the infrared dynamics of 
the theory. It is not necessary to include separate modes 
for infrared scales below Aqcd that occur in perturbation 
theory diagrams, as is shown explicitly for DIS. 

The role of the messenger scale in SCET will be 
clarified by studying a measurable quantity, the spin- 
dependent gi N structure function for DIS off a nucleon 
target. The structure function g\ has been chosen rather 
than the unpolarized structure functions because 
all the computations needed for the analysis of g\ have 
been given explicitly in the literature in the form needed 
here 111ill . The conclusions derived below are general, and 
apply equally well to other problems studied using SCET. 

The typical momentum scales in DIS at generic values 
of x are Q 2 and Aq C d . The final hadronic state X has 
invariant mass 

Ml = X —^Q 2 + M 2 (1) 

x 

where Mp is the target mass. As x —+ 1, M\ —* Q 2 ( 1 — 
x) < Q 2 , so the final hadronic state X becomes jet-like. 
DIS in the x —> 1 limit can be studied Hi using SCET. 
The large scale is Q , and the power-counting parameter A 
is taken to be 1 — x = A 2 <C 1. It is convenient to choose 
A 2 = Aqcd/Q, so that the jet scale is Q\ = -y/QAQCD- 
The standard factorization theorem for DIS i implies 
that 



where f/i q G are perturbatively calculable hard scatter¬ 
ing kernels, and /a<j/jv, Iag/n are the polarized quark 
and gluon distributions in the nucleon. The parton dis¬ 
tribution functions depends on the target, but the hard 
scattering kernels do not. As x —> 1, the hard scatter¬ 
ing kernels can be written as the convolution of a hard 
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FIG. 1: Box diagram contributing to deep-inelastic photon- 
gluon scattering via the process 7 * + g —> q + q. The total 
cross-section is given by the imaginary part of the diagram. 
The wavy lines are the high-Q 2 virtual photon, and the helical 
lines are the gluon target. 



FIG. 2: Graph contributing to the gluon matrix element of 
the quark distribution function. The two ® and the dashed 
line denote the insertion of the bilocal operator 0/\ q defined in 
Eq. 0. Graphs where the gluons connect to the Wilson line 
(the dashed line) vanish for the polarized parton distribution 
function. 


function with a jet function, and a soft function can be 
factored out of /a 9 ,ag QLlZl- This additional factoriza¬ 
tion will not be important in what follows. The stan¬ 
dard factorization form of Eq. <©> writes g\ as a sum of 
terms of the form K ® /, where K depends on Q but 
not Aqcd, and / depends on Aqcd but not Q. (® de¬ 
notes the convolution.) The main conclusion of Ref. Q is 
that the factorization form Eq. © needs to be modified 
so that the parton distribution functions also depend on 
Aqcd \J Aqcd / Q, and hence on Q, so that even though 
the convolution form Eq. <©> holds, the Q and Aqcd de¬ 
pendence has not been factorized. 

The hard scattering kernels are target independent, 
and so can be computed from a perturbative computa¬ 
tion of DIS off a partonic target. In particular, DIS off a 
gluon target can be written as 


9ig = 9lq ® fAq/G + 9ig ® /ag/g- (3) 


The graphs contributing to DIS scattering off a gluon 
target are the box graph in Fig. 0 and related graphs 
given by permuting the external lines. Evaluating the 
graphs explicitly in the Q —> 00 limit for an external 
gluon state with p 2 7 ^ 0 and a quark with mass m and 
unit electric charge gives 00 


91G 


4-7T 


( 2x — 1 ) In 


+3 — 4x + 


<3 2 (i - x) 


m 2 x — p 2 x 2 (l — x) 

p 2 x{ 1 — x) 
m 2 — p 2 x( 1 — x) 


(4) 


where h is the helicity of the gluon. This result shows 
that, in the limit x —> 1 , the graph depends on the scales 
Q 2 (1 — x) , m 2 , and p 2 (1 — x) . The mass to and p are taken 
to be infrared quantities of order Aqcd, so p 2 (l — x) is 
the messenger scale Aqcd \J Aqcd /Q in the limit x —> 1 . 

Eq. © can be used to compute the hard scattering 
kernel (Jig- In Eq. ©, (jiG- and fAq/G are °f order a s , 
so one only needs the order one expressions 


Substituting in Eq. m gives 

9ig = -^fAq/G + hgi G + O (a 2 ), (6) 

so to obtain g\ G requires a definition of fAq/G- 

The parton distributions are defined as the matrix el¬ 
ements of light-cone correlation functions of quark and 
gluon fields [jQj. The polarized quark distribution is the 
matrix element of the bilocal operator [r| 

C>Aq(k + ) = J dz~e~ lz k+ [ijj(z~)W(z~)'y + ^ 5 ijj(0) 

+4>(0)W\z~)j + 'y 5 ip(z~)\ (7) 

where ± denote light-cone directions, and W{z~) is a 
light-like Wilson line from 0 to the point y with light- 
cone coordinates y + = 0, yj_ = 0, y~ = z~. 

The parton distribution f Aq/T is ih e matrix element 
of OAq in the target T with momentum P, 

fAq/ T {x) = (T,P\0 Aq (xP + )\T,P). ( 8 ) 

Using an off-shell gluon target in Eq. m gives 0 (see 
Fig.© 


fAq/G ~ 


«s 

27T 


(2x — 1 ) In 
-1 + 


9 


m 2 — p 2 x{ 1 — x) 
2 


TO 


TO 2 — p 2 x( 1 — X) 


(9) 


which for to = 0 reduces to 


fAq/G = h 


27T 


( 2x — 1 ) In 


—p 2 x( 1 — x) 


- 1 


• (10) 


Substituting Eq. © in Eq. © gives 


3i G = 


4-7T 


(2x — 1 ) In 


Q 2 { 1 - x) 

p 2 x 


+ 3 — 4x 


• ( 11 ) 


fAG/G = hS( 1 - X), 


9U 


= 2^(1 ~ x )‘ 


The computation of the kernel g\ G is analogous to a 
(3) matching computation in an effective field theory. One 
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computes amplitudes in the full and effective theories, 
each of which is infrared sensitive, in perturbation theory, 
and takes the difference to obtain the values of matching 
coefficients, which are independent of the infrared scales. 
Here, the structure function g\ G and the parton distri¬ 
bution function f Aq /G are both sensitive to the infrared 
regulators m and p 1 but this dependence cancels in the 
kernel gi G . In the limit x —+ 1, 


9ig 


a s 

4-7T 


In 


Q 2 ( 1 - x) 



( 12 ) 


and depends on the jet scale Q 2 (l—x) ~ QAqcd, but not 
on any lower scales such as Aqcd or Aqcd \J Aqcd / Q , 
as expected. In SCET, the hard function is given by the 
matching coefficients at the scale Q, and the jet function 
is given by integrating out collinear modes at the scale 
Q 2 ( 1 — x). The Sudakov double-logarithms as x —> 1 
are summed by integrating SCET anomalous dimensions 
between Q 2 and Q 2 ( 1 — x) mum- All of these are 
perturbation theory computations, and are valid because 
the scales Q 2 and Q 2 (l — x) are both much larger than 
Aqcd. 

That the hard scattering kernel Eq. m is indepen¬ 
dent of infrared scales is usually taken as a proof of fac¬ 
torization. However, for factorization to hold, we also 
need the parton distributions to be independent of Q. 
An examination of the parton distribution Eq. II1UI1 for 
m = 0 shows that when 1 — x ~ Aqcd/Q, f Aq /G depends 
on the messenger scale Aqcd \J Aqcd / Q, and hence im¬ 
plicitly on Q , violating factorization. An evaluation of 
the diagram in Fig. fusing the method of regions ITol 
shows that the messenger region, with momentum of or¬ 
der p 2 (l — x ) contributes to the graph Jj, and it is this 
region that lead to the Q dependence in Eq. m- How¬ 
ever, the parton distribution Eq. © for to ^ 0 does not 
depend on the scale p 2 (l — x) for x —> 1, since in this 
limit to 2 p 2 (l — x) and the logarithm depends only on 
m 2 . 

It is important to remember that neither the compu¬ 
tation of gi G , nor that of f Aq /c is a valid computa¬ 
tion. Both depend on infrared physics for which per¬ 
turbation theory is not valid. Thus the dependence of 
gi G and f Aq /c on p 2 (l — x), and the necessity of messen¬ 
ger modes to reproduce Eqs. CTTU 1 ) is irrelevant. What 
has been determined so far is that the kernel is given 
by Eq. 11II) . and this can be obtained by a perturbative 
effective field theory computation using SCET. The de¬ 
pendence of the parton distribution on infrared scales is 
determined using Eq. m for a nucleon target. The rele¬ 
vant non-perturbative quantities it can depend on are the 
masses of physical hadrons that couple to the nucleon via 
the axial current matrix element. The spectrum of these 
hadrons cannot be computed in perturbation theory, but 
is known from both experimental observation and the¬ 
oretical considerations to not contain states arbitrarily 
close in mass to the nucleon. The nucleon parton dis¬ 
tribution depends on the light-cone Fourier transform of 


the nucleon matrix element of the gauge invariant bilocal 
operator {[){z~)W{z~)^ + 'y 5 ip(d). The z“ dependence of 
this correlation function is governed by the scale Aqcd 
in QCD, because of confinement, and so the parton dis¬ 
tribution only depends on this scale. Gauge invariance 
is crucial for the argument, since there can be massless 
gauge variant modes even in a confining gauge theory. 
That the operator is bilocal, rather than local, is not im¬ 
portant. Confinement implies there are no long range 
correlations even for non-local observables such as Wil¬ 
son loops. 

The above argument can be made more explicit by 
using the gauge invariant bilinears 

’J'O") = 4>i{z~)ip{z ~), (13) 

where </>„ is a color-triplet field with Lagrange density 
4>n(in • D)cj) n . The sole purpose of cj> n is to reproduce 
the Wilson line W(z~), so there is no longer an explicit 
Wilson line between 0 and z _ . In terms of these gauge 
invariant bilinears, 

0 Aq (k+) = ±-j dz-e~ iz ~ k+ [’F(z-) 7 + 75 ^( 0 ) 

+(0)7+75T(z - )]. (14) 

Inserting a complete set of states X between the d' 
and ik fields, which is now possible because there is no 
Wilson line, shows that 0 Aq {xP + ) couples the nucleon to 
gauge invariant intermediate states with Pjt = P+(l— x). 
These intermediate states are color singlet hadrons made 
of two quarks and one color-triplet </>„ field. The (t> n field 
is equivalent to a massless color-triplet quark moving in 
the light-cone + direction. The parton distribution func¬ 
tion depends on the spectrum of these excited states in 
QCD, which is governed by Aqcd- Thus the nucleon par- 
ton distribution function is not sensitive to scales below 
Aqcd- Note that there can be excited states with individ¬ 
ual mass differences which are accidently much smaller 
than Aqcd ■ In the DIS regime, one sums over many final 
states, and what matters is the typical mass splitting, of 
order Aqcd- In inclusive B decays, which are similar 
in structure to DIS, the relevant non-perturbative scale 
has been shown to be Aqcd 0 , even though there are 
individual states, such as the B*, which have an order 
Aqcd ! m b nrass difference from the B. 

It is instructive to see why factorization fails at the 
level of perturbation theory. The bilocal operator Eq. © 
is not Lorentz invariant, but is invariant under boosts 
along the z-axis. It depends on the momentum fc+, and 
a renormalization scale /x, but is independent of the large 
momentum q^. The matrix element of the bilocal op¬ 
erator can depend on the boost invariant combination 
p~p\ = p 2 { 1 — x). The perturbative computation of the 
parton distribution Eq. © depends on the spectrum of 
intermediate physical states in perturbation theory. Con¬ 
sider first the case of massless fermions, with p 2 ^ 0. This 
is the infrared regulator used in Ref. Q and in all the ex¬ 
amples where messenger modes contribute. The interme¬ 
diate states in Fig©are 4> n q states with massless quarks, 
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and can have any invariant mass > 0. There are interme¬ 
diate states with mass p 2 (l — x) = Aqcd \J Aqcd /Q, and 
so the parton distribution function depends on this scale, 
as can be seen from Eq. © with to 2 = 0. If to 2 ^ 0, then 
the lowest intermediate state has invariant mass to 2 . In 
the limit x —> 1 , to 2 — p 2 x(l — x) — > to 2 , and the par- 
ton distribution is sensitive to the scale to, but not to the 
lower scale p 2 (l—x). The existence of the messenger scale 
Aqcd \J Aqcd /Q depends on the precise form of the in¬ 
frared regulator Messenger modes are needed if one 
wants to reproduce the DIS cross-section for an off-shell 
partonic gluon target with massless quarks. However, we 
are interested in the DIS cross-section off hadrons, not 
partons; the partonic cross-section was merely an inter¬ 
mediate step in the computation of the kernel <fi G . The 
existence of low-mass hadronic states is an artefact of 
perturbation theory. Low-mass states are not present in 
QCD, messenger modes are not needed, and factorization 
holds. 

There is an alternate argument for messenger modes 
in DIS as a; —> 1. The computation of DIS in the 
Breit frame shows that there are contributions from h- 
collinear modes, n-collinear modes and ultrasoft modes. 
The n-collinear mode describe the target, the n-collinear 
mode the final hadronic jet, and the ultrasoft modes 
the low energy gluons. For the matching at the hard 
scale Q and the running between Q and the jet scale 
Q 2 (l — x) ~ Q Aqcd, one can treat both the n and h 
collinear modes as having an off-shellness (for the pur¬ 
poses of regulating the infrared divergences and power 
counting) QAqcd- In this case, the momentum of these 
modes in the Breit frame is ~ Aqcd, Pn ~ Q, 
Pn ~ Qi Pn ~ Aqcd, and the ultrasoft gluons have 
typical momentum p 2 ~ p+Pn ~ Aq C d . The target 
rest frame is obtained from the Breit frame by an or¬ 
der Q/ Aqcd boost along the collision axis. If, in addi¬ 
tion, one lowers the off-shellness of the n-collinear mode 
Aq CD @, the ultrasoft modes in the Breit frame get 
mapped into messenger modes in the target rest frame 
with p 2 ~ Aq CD /Q. These messenger modes are the ana¬ 
log of rescaled ultrasoft modes in the Breit frame which 
also have p 2 ~ Aq CD /Q . However, an explicit compu¬ 
tation shows that this is not the case jjjj. The ultrasoft 
modes in the target rest frame reproduce the contribu¬ 
tions of both the n-collinear and messenger modes in the 
Breit frame. Equivalently, the parton distribution func¬ 
tion in Eq. 0 involves only ultrasoft fields, and a di¬ 
rect computation of its matrix element 0 gives Eq. ©. 
Fig© gets contributions from regions where some quark 
lines have invariant mass of order p 2 , and others have 
invariant mass of order p 2 ( 1 — x)] however, both are in¬ 
cluded in the ultrasoft graph. 

In the target frame computation, the target quarks 
have an off-shellness of order Aq CD from the very be¬ 
ginning, and are described by ultrasoft fields. Once the 
jet scale has been integrated out, the dynamics of these 
ultrasoft fields gives the non-perturbative parton distri¬ 
bution function. In the Breit frame computation, the 


n-collinear and ultrasoft/messenger modes are both non- 
perturbative below the jet scale, and together reproduce 
the effects of the non-perturbative ultrasoft fields in the 
target rest frame. The non-perturbative effects of both 
Breit frame modes can be combined into a single non- 
perturbative function, which has no knowledge of the 
messenger scale. Modes such as the messenger mode arise 
if one tries to “over-factor” a cross-section, i.e. separate 
out all the different scales that appear in a perturbation 
theory graph into separate contributions, even if all the 
scales are non-perturbative. In DIS, they would appear if 
one tries to factor the non-perturbative parton distribu¬ 
tion into separate non-perturbative pieces based on the 
momentum scaling of perturbation theory. 

Note that Ref. [§] did not attempt to compute the to¬ 
tal cross-section for scattering off a partonic target; par¬ 
tonic matrix elements were used only to compute match¬ 
ing coefficients and anomalous dimensions. In these 
computations, for which perturbation theory is valid for 
all the dynamical modes, one can treat the n-collinear 
and ultrasoft modes separately in the Breit frame. The 
matching conditions and anomalous dimensions are iden¬ 
tical in the Breit and target rest frame. The two Breit 
frame modes below the jet scale, where they become non- 
perturbative, were combined together in a single non- 
perturbative function. 

It might be possible to factor a soft function out of 
the parton distribution / = S (g> /, as suggested by the 
Breit frame calculation, where S is the matrix element of 
Wilson line operators, and / is the matrix element of n- 
collinear fields (see e.g. Ref. Q). The parton distribution 
/ renormalized at a low scale p is independent of Q, 
and any definition of the split of / into S and / in the 
theory below the jet scale, to be useful, should also be 
independent of Q , and so independent of the messenger 
scale. 

SCET is an effective theory formulated with ultrasoft 
and collinear modes H, H, H, IH , and provides an adequate 
description of DIS as x —> 1 off hadronic targets. It al¬ 
lows one to separate short distance contributions at the 
hard scale Q 2 and the jet scale Q 2 ( 1 — x) from the non- 
perturbative contributions. It is not necessary to have 
separate propagating modes for every momentum scale 
that contributes to a perturbative diagram with some in¬ 
frared regulator. In the target rest frame, for example, 
all non-perturbative effects are included in the ultrasoft 
contribution, and messenger degrees of freedom are not 
needed. In general, one does not need a propagating 
mode in an effective field theory for every momentum 
scale in the problem with some infrared regulator. A 
simple example is NRQED, with expansion parameter 
a. The Hydrogen and Positronium energy levels can be 
computed to high order using NRQED, which has poten¬ 
tial {E ~ mct 2 ,p ~ too), soft (E ~ majp ~ met) and 
ultrasoft (E ~ ma 2 ,p ~ met 2 ) modes [20|. The shape 
function for Positronium decay, which is analogous to the 
structure function in DIS, has been computed explicitly 
using this theory {^]. It is not necessary to introduce 
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additional modes at the scale ma 4 of the hyperfine split¬ 
ting, even though the ortho-para Positronium mass split¬ 
ting Am = 7m e a 4 /12 occurs in energy denominators of 
the form Am — E 1 . 
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